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Abstract 

The ATLAS Collaboration has measured the inclusive production of Z bosons via their de¬ 
cays into electron and muon pairs in p +Pb collisions at \Js NK = 5.02 TeV at the Large 
Hadron Collider. The measurements arc made using data corresponding to integrated lumi¬ 
nosities of 29.4 nb -1 and 28.1 nb -1 for Z —» ee and Z —> pp, respectively. The results from 
the two channels are consistent and combined to obtain a cross section times the Z —> ££ 
branching ratio, integrated over the rapidity region \y* z \ < 3.5, of 139.8 ± 4.8 (stat.) ± 6.2 
(syst.) ± 3.8 (lumi.) nb. Differential cross sections are presented as functions of the Z 
boson rapidity and transverse momentum, and compared with models based on parton dis¬ 
tributions both with and without nuclear corrections. The centrality dependence of Z boson 
production in p +Pb collisions is measured and analyzed within the framework of a standard 
Glauber model and the model’s extension for fluctuations of the underlying nucleon-nucleon 
scattering cross section. 


© 2015 CERN for the benefit of the ATLAS Collaboration. 

Reproduction of this article or parts of it is allowed as specified in the CC-BY-3.0 license. 




Contents 


1 Introduction 2 

2 The ATLAS detector 3 


3 Analysis 4 


3.1 

Data sample 

4 

3.2 

Lepton reconstruction 

4 

3.3 

Centrality 

5 

3.4 

Monte Carlo simulation corrections 

7 

3.5 

Yield extraction 

8 

3.6 

Systematic uncertainties 

10 

3.7 

Lepton channel combination 

11 

Results 

12 

4.1 

Z —> U cross section 

12 

4.2 

Centrality-dependent yield 

14 

Summary 

18 


1 Introduction 


The study of electroweak bosons in Pb+Pb collisions at the Large Hadron Collider (LHC) at CERN 
has demonstrated that the production rate of non-strongly interacting particles scales with the number of 
nucleon-nucleon collisions, A co ii. This has been observed for photons [1], W bosons [2, 3], and Z bosons 
[4, 5]. The momentum and rapidity distributions of Z boson yields are consistent with Pythia [6] simu¬ 
lations of pp collisions multiplied by the average nuclear thickness function, (Taa)> which is equivalent 
to (iV C oii> divided by the total nucleon-nucleon cross section [4]. Z boson production in Pb+Pb collisions 
was found to be consistent with Next to Leading Order perturbative Quantum Chromodynamics (NLO 
QCD) calculations that disregard nuclear modifications in the treatment of parton distribution functions 
(PDF). However, nuclear modification is not excluded within the precision of the measurement [7]. The 
production of Z bosons, when examined as a function of centrality, was also found to scale with (/V co |]>. 

To differentiate between initial- and final-state effects in heavy ion (HI) collisions, the study of p +Pb 
collisions is used at the LHC. One could expect that the hot and dense QCD medium cannot be formed 
in such collisions, unlike in the Pb+Pb case, and that modifications to the final-state particles relative to 
nucleon-nucleon collisions should originate from the initial state of the nucleus. This assumption was 
challenged by the very first results from p +Pb collisions at ^Js NN = 5.02 TeV produced at the LHC 
in 2012. Results on multi-particle correlations, published by three LHC experiments [8-14], revealed 
collective behavior in p +Pb collisions similar to that previously measured in HI collision systems. The 
yields of jets measured by ATLAS scale with (A co ii) when measured inclusively for all centralities, but 
show significant deviations from binary scaling when considered in centrality selections [15]. The CMS 
Collaboration has measured dijet pseudorapidity distributions and found them to agree better with pre¬ 
dictions that include nuclear PDF modifications than with predictions that do not include nuclear effects 
[16]. The CMS Collaboration has also recently measured production of W bosons in p +Pb collisions and 
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observed hints of nuclear modifications of the PDF [17]. Collectively, these results have highlighted the 
need for a better understanding of the initial conditions of p +Pb collisions. 

Unlike symmetric Pb+Pb collisions, in p +Pb collisions nuclear modifications of the PDFs in the lead nu¬ 
cleus create an asymmetry in the rapidity-dependent cross section of Z bosons; this presents an attractive 
observable for the study of initial-state nuclear conditions. The centrality-dependent yield of Z bosons is 
a well suited probe to test our understanding of p +Pb collision geometry. The LHCb Collaboration has 
made a first exploratory measurement of Z bosons at far forward and backward rapidities [18] based on 
an integrated luminosity of 1.6 nb _1 . 

This paper presents the results of the measurement of Z boson production in p +Pb collisions at y.s\ vv = 
5.02 TeV using the ATLAS detector. The yield of Z bosons is measured as a function of their transverse 
momentum pj, rapidity in the center-of-mass frame (y* 7 ), 1 and centrality. The leptonic decays of the 
Z boson (Z —> ee and Z —> pp) are used for its reconstruction. In the muon channel it is possible to 
reconstruct the Z boson in the rapidity range -3 < if z < 2, while in the electron channel this range can be 
extended to \y * z | <3.5. The larger acceptance in rapidity for Z —> ee candidates is possible because of the 
larger acceptance of the ATLAS calorimeters compared to the muon spectrometer (MS; Secs. 2 and 3.2). 
The efficiency of Z boson reconstruction is calculated from detector simulations (Sec. 3.4). Backgrounds 
in each channel are estimated using simulations and data-driven methods (Sec. 3.5). Results measured 
in the dimuon and dielectron decay channels are combined after accounting for uncertainties and their 
correlations (Sec. 3.7). The measured cross sections and centrality-dependent yields are compared to 
models of p +Pb collisions composed of 82 pp and 126 pn collisions in which the production of Z bosons 
is obtained using perturbative QCD calculations. 


2 The ATLAS detector 

The ATLAS detector [19] at the LF1C covers nearly the entire solid angle around the collision point. It 
consists of an inner hacking detector surrounded by a thin superconducting solenoid, electromagnetic and 
hadronic calorimeters, and the MS. 

The inner-detector (ID) system is immersed in a 2 T axial magnetic held and provides charged particle 
tracking in the pseudorapidity range I 77 I < 2.5. It comprises a high-granularity silicon pixel detector cov¬ 
ering the collision region, surrounded by a silicon microstrip tracker and a transition radiation tracker. 

The calorimeters cover the range |//| < 4.9. Within the region |//| < 3.2, electromagnetic calorimetry is 
provided by barrel and endcap high-granularity lead liquid-argon (LAr) calorimeters, with an additional 
thin LAr presampler covering \i]\ <1.8. Behind the electromagnetic calorimeter there is a steel/scintillator 
sampling hadronic calorimeter covering |//| < 1.7, and LAr hadronic calorimeters extend the coverage to 
|//| < 4.9. Forward electromagnetic calorimeters (FCals) are located in the range 3.1 < |?/| < 4.9. Electrons 
may be reconstructed over the entire electromagnetic calorimeter system, |//| < 4.9. 


1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the center of the LHC ring, and the i/-axis points upward. 
Cylindrical coordinates (r, <f>) are used in the transverse plane, (j> being the azimuthal angle around the beam pipe. The rapidity 
in the laboratory frame is given by y ub = | In and pseudorapidity is defined as 77 = - ln[tan(0/2)]. Positive rapidity 
corresponds to the direction the proton beam travels, “proton-going", and negative rapidity is referred to as “Pb-going". In 
this convention the asymmetric beam energies result in a center-of-mass shifted to rapidity y* = y klb - 0.465. 
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The MS comprises separate trigger and high-precision tracking chambers that measure the deflection 
of muons in a magnetic field generated by superconducting air-core toroids. The precision chambers 
cover the region |//| < 2.7 with three layers of monitored drift tubes (MDT) chambers, complemented by 
cathode strip chambers (CSC) in the innermost layer of the forward region. The muon trigger system 
covers the range |p| < 2.4 with resistive plate chambers (RPC) in the barrel (It/I < 1.05), and thin gap 
chambers (TGC) in the endcap regions (1.05 < M < 2.4). 

The ATLAS detector has a three-level trigger system [20]: the hardware-based level-1 (LI) trigger and the 
software-based High Level Trigger (HLT), which is subdivided into the Level-2 (L2) trigger and Event 
Filter (EF). Single-electron and single-muon triggers are used to acquire the data analyzed in this paper. 
Minimum-bias events are selected based on signals in the minimum-bias trigger scintillators (MBTS) that 
detect charged particles in the range 2.1 < \r/\ < 3.9. 


3 Analysis 

3.1 Data sample 

This analysis uses the 2013 ATLAS p +Pb collision data at ^/s NN = 5.02 TeV, produced from a 4 TeV 
proton beam and a 1.57 TeV per nucleon lead beam. The asymmetric energy of the beams resulted in a 
shift of the center-of-mass by 0.465 units of rapidity relative to the laboratory frame. After 60% of the 
data were recorded the directions of the proton and lead beams were reversed. Results obtained in the 
two data periods are found to be consistent with each other. In this paper all data from both periods arc 
presented using the convention that the proton beam travels forward in the positive rapidity direction. 

Following stringent data-quality requirements, the Z —> ee and Z —> pp analyses use samples correspond¬ 
ing to integrated luminosity values of 29.4 ± 0.8 nff 1 and 28.1 ± 0.8 nbr 1 , respectively. The luminosity 
measurement for the 2013 p+Pb data is calibrated based on dedicated beam-separation scans. System¬ 
atic uncertainties similar to those studied for the calculation of pp luminosity [21] are calculated. The 
combination of these systematic uncertainties results in a total uncertainty in the ATLAS luminosity scale 
during proton-lead collisions at *\Js NN = 5.02 TeV of 2.7%. Minimum-bias p+Pb collisions were selected 
by a trigger based on a signal in both MBTS counters. Minimum-bias events are required to have the time 
measured in each MBTS be consistent within 10 ns, and a reconstructed collision vertex within 175 mm 
of the nominal collision point in the longitudinal direction [ 22 ]. 


3.2 Lepton reconstruction 

Electron candidates are first identified by the LI trigger as a cluster of cells in the electromagnetic 
calorimeter, formed into (A <p x A 77 ) = 0.1 x 0.1 trigger towers, within the range |//| < 2.5 and with 
the cluster transverse energy exceeding 5 GeV. The HLT then incorporates tracking information from the 
ID and imposes electron identification requirements on the electron candidates. A trigger for at least one 
electron candidate with Ej > 15 GeV and satisfying loose identification requirements is used to select 
events. 

In the offline analysis, electron candidates within |p| < 4.9 are selected using the ATLAS reconstruction 
algorithm [23]. Electrons with |p| < 2.47, referred to as midrapidity electrons, require the matching of 
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a track to an energy cluster in the electromagnetic calorimeter. In addition to the reconstruction require¬ 
ments, further electron identification selections based primarily the shower shape in the electromagnetic 
calorimeter are made to reject background. Electron identification requirements used in previous ATLAS 
analyses [23] are used to provide quality classification of electrons based on the tightness of the identi¬ 
fication criteria they satisfy. The triggering electron is required to have Ej > 20 GeV, to be outside the 
pseudorapidity interval 1.37 < |//| < 1.52 (a transition region between the barrel and endcap calorimeters 
which contains a relatively large amount of inactive material), and to satisfy tight identification quality re¬ 
quirements. If the other electron is within |//| < 2.47, it must have Ej > 10 GeV and satisfy looser quality 
requirements. Forward-electrons are those reconstructed within the range 2.5 < |//| < 4.9 based on energy 
deposited in the FCal [23]. There is no tracking in this region, so the electron candidate reconstruction 
and identification is derived solely from the calorimeter signal and does not have an associated charge. 
Forward-electrons are required to have Ej > 20 GeV. Of Z ee decays with \if 7 \ < 3.5, approximately 
82% fall into the fiducial acceptance defined by the electron pj and t) requirements. 

Muon candidates are first identified at the LI trigger, based on hits in either the RPC or TGC. The HLT 
then reconstructs muon tracks in the vicinity of the detector region reported by the LI trigger. The L2 
trigger uses an algorithm to perform a fast reconstruction of muons, which is then refined in the EF by 
incorporating the hits in the ID tracking as well as those in the MS tracking. Events containing at least 
one muon with p-\ greater than 8 GeV are accepted by the HLT. 

For the Z —> pp analysis, muons are identified from candidates reconstructed in both the MS and ID [24]. 
Muons are reconstructed separately in the MS and ID and a ^-mini miz ation procedure is used to obtain 
combined muon kinematic information. To reduce background from jets, each muon is required to pass a 
loose track-based isolation selection. Tracks are considered in a cone of size A R = \/(A//) 2 + (A<f>) 2 = 0.3 
around the direction of the muon. The muon is considered isolated if the scalar sum of the p-\ of these 
tracks, excluding the muon, is less than 50% of the muon pj. The efficiency of this selection is greater 
than 99%. The triggering muon is required to have pj > 20 GeV and be within |;/| < 2.4. The second muon 
must be within \p\ < 2.47 and have pj > 10 GeV. Approximately 68% of Z bosons with -3 < if z < 2 fall 
into the fiducial acceptance defined by the muon kinematic requirements. 


3.3 Centrality 

In addition to measuring the Z boson cross section, the Z boson yield per minimum-bias event is measured 
for different centrality selections. In order to characterize the p +Pb collision geometry, each event is 
assigned a centrality based on the total transverse energy measured in the FCal on the Pb-going side 
(-4.9 < rj < -3.2), Z£' , | ' CaI [22]. Collisions with more (fewer) participating nucleons are referred to as 
central (peripheral). As in Ref. [22], the standard Glauber model [25] approach is used to calculate the 
mean number of participating nucleons, (A pa rt). The mean number of inelastic nucleon-nucleon collisions, 
(Vcoii), is (Apart) _ 1- Based on the observed centrality dependence of the charged particle multiplicity 
[22], a Glauber-Gribov Color Fluctuation (GGCF) model [26, 27], an extension to the Glauber model 
which allows event-by-event fluctuations of the nucleon-nucleon cross section cr(N + N —> X), is also 
considered. In this model the magnitude of the fluctuations is characterized by the parameter oj (r , with 
o) ir =0 corresponding to the standard Glauber model. Following Ref. [22], two values of <o (T , 0.11 and 0.2, 
based on the calculations in Refs. [26, 27], are implemented and considered for this analysis. 

Following Ref. [22], the centrality selections used for this analysis, in order from most central to most 
peripheral, are 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-60%, and 60-90%. For the study of the 
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y* z distribution as a function of centrality (see Sec. 4.2), larger bins are used: 0-10%, 10-40%, and 40- 
90%. For the most peripheral collisions, centrality greater than 90%, centrality modeling, and associated 
geometric quantities are not well constrained. Pileup events, those containing multiple p +Pb interactions 
from the same bunch-crossing, are removed by rejecting events in which more than one primary vertex 
is reconstructed. The fraction of candidate events which is removed from the centrality-selected Z boson 
yield analysis due to pileup rejection is approximately 5%. Diffractive events are identified by a rapidity 
gap (defined by the absence of calorimeter energy clusters) of greater than two units on the Pb-going side 
of the detector, and excluded. This leads to a rejection of less than 0.1% of Z boson candidate events. 
The total number of minimum-bias events corresponding to the luminosity sampled by the trigger, /V cvt , 
is used to define the Z boson yield per event in each centrality selection. 

Besides its sensitivity to the event geometry, the X£^ Cal for a fixed geometry may also be affected by 
the presence of a hard scattering process in the event. In particular, the calculation of centrality for 
Z —> ee events in which there is a forward-electron in the Pb-going side FCal, is biased by the energy 
of the electron. This is corrected by subtracting the transverse energy of the electron from EE£ Cal . The 
subtraction procedure is found to effectively recover the correct centrality of minimum-bias events into 
which simulated Z —> ee decays containing electrons in the Pb-going side FCal were overlaid. 

In addition to the case where there is a Z-decay electron in the Pb-going side FCal, Z/f , | ' < al may be more 
subtly biased in all Z boson events. The presence of a Z boson (or any hard process) is correlated with 
a higher transverse energy of the underlying event. Consequently, more energy may be deposited in the 
Pb-going side FCal in events containing a hard scattering process than in those that do not contain one. 
This causes a bias in the centrality-dependent yield, as the Z boson yield is enhanced in the more central 
events but depleted in the more peripheral ones. This effect, referred to as a centrality bias, has been noted 
for yields of hard processes in <r/+Au collisions at yjs NN = 200 GeV by the PHENIX Collaboration [28], 
and by the ALICE Collaboration in p +Pb collisions at ^Js NN = 5.02 TeV by the ALICE Collaboration 
[29], 

A correction to the centrality-dependent yields of hard processes in c/+Au collisions at tJs nn = 200 GeV 
has been studied by the PHENIX Collaboration [28] (the correction is also calculated for p +Pb at -\]s NN = 
5.02 TeV). The centrality bias correction used by the PHENIX Collaboration is based on the modeling of 
an increase in the mean particle multiplicity produced by the specific NN collision that undergoes a hard 
scattering. Recently, similar calculations of a centrality bias have been made in which all NN collisions 
may contribute to an increase in the particle multiplicity [30]. The increase in multiplicity stemming from 
each NN collision is taken to be proportional to the contribution from that collision to the Ej in the event. 
This model is applied to the ATLAS p +Pb centrality classification for the standard Glauber analysis as 
well as the GGCF models, and thus used to calculate corrections to the hard process yield measured in 
a given centrality bin. Because the /V par[ probability distribution varies less steeply in the GGCF models 
than in the standard Glauber model, the centrality bias corrections are closer to unity for the GGCF cases. 
The corrections from Ref. [30] are shown in Figure 1 . The reciprocals of the corrections are applied as 
multiplicative factors to the centrality-dependent Z boson yields measured in the present analysis. 

Using Z bosons measured in ATLAS, data-driven centrality bias corrections may be calculated and com¬ 
pared with the results from Ref. [30]. To do so, the corrections are calculated by comparing the transverse 
energy deposited in the FCal in events selected by the minimum-bias trigger and in Z boson events from 
pp collisions (in which there is no centrality to consider). This effect is studied in pp collisions at 
-\fs = 2.76 TeV and \[\ = 7 TeV from the 2011 LHC run. At both energies a significant increase in 
the mean transverse energy deposited in the FCal is observed, and within the uncertainties of the mea¬ 
surement found to be independent of the Z boson kinematics. A single value is interpolated from the 
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Figure 1: (Color online) Centrality bias correction factors and their uncertainties (bars), from Ref. [30], for the 
Glauber and GGCF model configurations, and the bias correction factors derived from data (points) as explained in 
the text. The reciprocals of the correction factors are applied as multiplicative factors to the centrality-dependent Z 
boson yields. 


-^s = 2.76 TeV and yfs — 1 TeV data for -\js NN = 5.02 TeV. The interpolation is performed using 
both logarithmic and linear expressions, and the difference between them contributes to the systematic 
uncertainty. A correction is made to account for the shifted center-of-mass in the /;+Pb system which 
changes the effective FCal acceptance in r/ compared to pp. From this procedure an additive shift to 
S£'^ Cal of 2.0 + 0.5 GeV is calculated. In each p +Pb Z —a (t event this value is subtracted from Z£ , ^ Cal 
and the resulting value used to determine a corrected centrality of the event. The centrality-dependent 
yield may be constructed, according to the method described in Sec. 3.5, using both the subtracted and 
unsubtracted E£^ Cal values and their ratio is then comparable to a centrality bias correction factor that 
may be compared to those calculated in Ref. [30]. As shown in Figure 1, within the uncertainties this 
data-based method is compatible with the model calculations. 


3.4 Monte Carlo simulation corrections 

The trigger, reconstruction, and identification efficiencies of electrons and muons as well as the muon iso¬ 
lation efficiency are evaluated by a Monte Carlo (MC) simulation complemented by data-driven estimates 
of these quantities. Using the Powheg generator [31] (with the CT10 PDF [32]) interfaced to Pythia8 
[6] for simulation of the parton shower, approximately 10 million Z —> ee and 4 million Z —> pp events 
were simulated. The Z bosons were generated from pp and pn collisions, which were added together with 
weights 82/208 and 126/208, respectively, corresponding to the numbers of protons and neutrons in the Pb 
ion. The response of the ATLAS detector to the generated particles was modeled using GEANT4 [33, 34], 
Due to the dependence of electron identification and reconstruction efficiency on detector occupancy, the 
simulated Z —» ee events were overlaid with data events selected with the minimum-bias trigger, and then 
reconstructed. 

To cross-check the efficiencies calculated in the MC simulation, a ‘tag-and-probe’ technique is employed. 
A ‘tag’ is defined as a fully reconstructed high-quality triggered lepton, whereas the ‘probe’ is a lepton 
candidate to which triggering, reconstruction, or quality requirements are not applied. Using tag-and- 
probe pairs with an invariant mass met consistent with selection of Z bosons, the efficiency of the probe 
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with additional requirements is calculated. The mass window used depends on the background present 
in the probe sample, and ranges from 80 < mu < 100 GeV to 87 < mu < 95 GeV. For example, the 
electron trigger efficiency is measured from high-quality reconstructed electron probes selected without 
an a priori trigger requirement, and the MS reconstruction efficiency is measured from charged-particle 
tracks in the ID without an a priori MS signal requirement. The MC simulation is scaled to match the 
efficiencies determined with the data-driven tag-and-probe method. The factors used to scale the MC 
electron response are derived from the 2013 p+ Pb data set. Muon reconstruction is insensitive to the 
differences between 2013 p+ Pb and 2012 pp conditions and therefore the scale factors for muons are 
taken from Z —» pp events collected in the 2012 pp dataset [24]. The scale factors for the Z —> ee 
(Z —» pp) MC events deviate from unity by less than 5% (1%). 

The trigger efficiencies of electrons and muons depend on ?/, and have average values of approximately 
95% and 83%, respectively. The reconstruction and identification efficiency of the more stringently se¬ 
lected electrons is approximately 75%, whereas the reconstruction and identification efficiency for the 
looser-quality midrapidity electrons is approximately 88% and the forward-electron efficiency is about 
65%. These values depend on // and pj. The muon reconstruction efficiency is approximately 95%, 
depending on rj. 

Correction factors for the yields of Z ee and Z —> pp candidates are calculated from the MC simu¬ 
lations as functions of y* z , and p+ Pb centrality. These corrections take into account the cumulative 
losses due to trigger, reconstruction, and identification efficiency as well as the kinematic acceptance of 
the decay leptons. The correction is defined relative to all generated Z bosons within the mass window 
66 < mu <116 GeV. The total efficiency for reconstructing a produced Z boson, including acceptance, 
is approximately 55% for Z —> ee and 65% for Z —* pp. Following the subtraction of background (see 
Sec 3.5) and application of the correction factor, a corrected yield of Z bosons is obtained in each bin 
of Pj and if 7 . The uncertainty on the correction factor follows from the uncertainty of the data-driven 
tag-and-probe checks of the MC, primarily due to the relatively low number of tag-and-probe events in 
the data. The uncertainty associated with the lepton identification efficiency is the dominant uncertainty 
for both the Z —> ee and Z —> pp analyses. This uncertainty is approximately 10% for Z —> ee (rising 
as high as 15% for pairs including a forward-electron), and 1.5% for Z —> pp. The sizable uncertainty 
for Z —» ee is primarily driven by the limited size of the tag-and-probe 2013 p+ Pb dataset. The other 
uncertainties are typically less than 2% in both channels. 


3.5 Yield extraction 

To form Z —» ee candidates, all electrons found in triggered events are paired with each other. When both 
electrons are at midrapidity (|^| < 2.47), the unlike-sign charged pairs with an invariant mass satisfying 
66 < m ee < 116 GeV are accepted as signal Z boson candidates. The like-sign pairs in this window are 
used to estimate the combinatorial background, created primarily by jets. In total, 1647 unlike-sign pairs 
and 52 like-sign pairs are reconstructed. The like-sign pairs are composed of combinatorial background 
and Z —> ee decays in which one of the electrons has misreconstructed charge. The contribution from 
pairs with misreconstructed charge is estimated, using the MC simulation, to be half of the like-sign pairs, 
and the remainder is taken as an estimate of the background. Pairs made of one midrapidity electron 
and one non-triggering forward-electron have a larger contribution from background and so an invariant 
mass window of 80 < m ee < 100 GeV is used to select Z boson candidates. To facilitate combination 
of all Z —> ee candidates, an acceptance correction is made to account for the smaller mass window. 
No charge requirement is made for these candidates because the non-triggering electron is outside the 


acceptance of the ID and therefore does not have a reconstructed charge. There are 264 such candidates, 
of which an estimated 5% are background based on a fit of the invariant mass distribution. The fit is 
performed in the range 60 < m ee <120 GeV using a signal shape from the MC simulation, and several 
background parametrizations assuming exponential or polynomial descriptions of the background. The 
mass distributions of Z —> ee candidates are shown in Figure 2(a) and 2(b), along with the reconstructed 
MC simulation of the same quantity. The estimated background is subtracted from the signal candidates 
differentially in rapidity, transverse momentum, and centrality. 



60 70 80 90 100 110 120 

rn„ [GeV] 

Figure 2: (Color online) The dilepton invariant mass distributions in data and MC simulation, (a) Z —» ee candidates 
with both electrons at midrapidity {\rj\ < 2.47); (b) candidates in which the non-triggering electron is a forward- 
electron (2.5 < \rj\ < 4.9); (c) Z —> /./// candidates. The midrapidity Z —» ee and Z —» /jfi distributions are 
overlaid with the mass distributions of like-sign lepton pairs. For candidates in which the non-triggering electron 
is a forward-electron (2.5 < \r]\ < 4.9) the background, estimated based on a decomposition of the invariant mass 
distribution (see text), is shown. In all plots the vertical lines indicate the mass window within which candidates 
are defined, and the MC simulation is normalized to the data inside this region. 


A similar procedure is also followed to select Z —» q/i candidates with an invariant mass of 66 < m /(/ , < 
116 GeV. This selection yields 2032 unlike-sign charged candidates and 4 like-sign pairs; their mass 
distribution is shown in Figure 2(c). The MC simulation describes the data well in both lepton channels. 
The slight shift of the mass peak visible between the data and the simulation for dielectron events has 
only a very small effect on the calculation of corrections based on the MC simulation and is incorporated 
into the systematic uncertainty associated with electron reconstruction. 

Based on the like-sign pairs and MC simulation of charge misreconstruction, the uncertainty from the 
background subtraction is approximately 1 % in the Z —» ee channel for pairs in which both electrons are 
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at midrapidity; in pairs involving a forward-electron the uncertainty ranges from 5% to 20% based on fits 
of the invariant mass distribution. The background uncertainty in the Z —> pp yield is negligible. The 
largest source of correlated background in both lepton decay channels is the decay of Z —» rr events into 
dielectron and dimuon pairs. These are simulated and reconstructed just as Z —> ee and Z —* pp are, but 
are found to have a negligible contribution following the analysis procedures. 


3.6 Systematic uncertainties 

The dominant source of uncertainty in the Z —> ee measurement stems from imperfect knowledge of the 
efficiency of the electron identification requirements. The uncertainty is driven by the limited number of 
events available for the tag-and-probe analysis which had to rely on 2013 p +Pb collision data because 
the electron reconstruction performance changed due to detector conditions and occupancy compared to 
earlier pp collision data. The uncertainty is larger in pairs involving a forward-electron, and the sample 
has a lower purity than the sample of midrapidity electrons. Other electron uncertainties are significantly 
smaller and are associated with trigger efficiency, electron reconstruction efficiency and energy resolution, 
background subtraction (which becomes significant for forward-electrons), and charge misreconstruction. 
In addition a small uncertainty stems from possible differences between the simulated y* z distribution and 
the one measured in data. The Z —> ee systematic uncertainties depend on p|, y* z , and /;+Pb centrality 
and are summarized in Table 1 . 


Source 

Uncertainty Range [%] 

Electron ID 

6-14 

Electron Reconstruction 

1-3 

Electron Trigger 

1-2 

Background 

1-3 

MC y* z Shape 

0-2 

Forward-Electron Reconstruction 

4-15 

Forward-Electron Background 

2-10 


Table 1: Relative systematic uncertainties, in percent, associated with the measurement of Z —* ee. The uncertain¬ 
ties typically increase at the more forward rapidities. Background includes charge misreconstruction, and electron 
reconstruction includes resolution. The last two rows refer only to pairs where one of the electrons was recon¬ 
structed in the range 3.1 < \rj\ < 4.9. 


The conditions of muon reconstruction in the 2013 p+Pb collision data closely resemble those in pp 
collisions described in Ref. [24], The small uncertainties from the more abundant pp data are used in this 
analysis. An uncertainty of 1%, based on the performance of the muon reconstruction in high-pileup pp 
collisions, is associated with the scale factors to account for possible differences between the datasets. 
The uncertainties depend on p| and y* z . Table 2 summarizes the systematic uncertainties on the Z —> pp 
measurement. 

In addition to the Z boson measurement uncertainties, a 2.7% uncertainty is associated with the luminosity 
calculation. For the centrality-dependent yields that arc scaled by (7\a) the uncertainties of the Glauber 
model calculations are taken from Ref. [22]. The correction to centrality due to bias from the presence of 
a hard scattering, taken from Ref. [30], has uncertainties as shown in Figure 1. 
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Source 

Uncertainty Range [%] 

Muon ID & Reconstruction 

1.5 

Muon Trigger 

1-2 

Background 

<1 

Pt Resolution 

<1 

MC y* z Shape 

<1.5 


Table 2: Relative systematic uncertainties, in percent, associated with the measurement of Z —> pp. 


3.7 Lepton channel combination 

The cross section in each leptonic decay channel is defined for the mass window 66 < my < 116 GeV, 
the rapidity ranges \y* z \ <3.5 for Z —» ee and -3 < y* z < 2 for Z —» pp, and the full decay-lepton 
kinematic phase-space. The Z —> ee and Z —> pp yields, corrected for acceptance and efficiency, are 
used to calculate the cross section in each channel, and good agreement between the two is observed as 
is shown for the y* z distributions in Figure 3. 



Figure 3: (Color online) (a) Differential Z boson production cross section, dcr/dy* z , as a function of Z boson rapidity 
in the center-of-mass frame y* z , for Z —* ee and Z —» pp (upper panel), (b) Their ratio. Bars indicate statistical 
uncertainty and the shaded boxes systematic uncertainty. 


The two decay channel results are combined to one set of Z —» £( data using the method described in 
Refs. [35, 36]. The technique uses a^ 2 minimization procedure with a nuisance parameter formalism 
to combine the data sets coherently. The procedure distinguishes between those systematic uncertainty 
sources that are uncorrelated bin-to-bin, uncorrelated across data sets, and fully correlated bin-to-bin and 
across data sets. In this way, combined points are calculated to optimize the overall agreement of the data 
sets, given the correlation of the uncertainties. This may result in differences in the combined Z —> tl 
data points relative to the Z —> ee data points in the rapidity regions in which there are no Z —» pp data 
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points. Following this, an integrated cross section for the region \y* z \ < 3.5 is defined for the combined 
Z —» It points based on both the Z —> ee and Z —* pp data even though the Z —> pp data arc limited to 
-3 < y* z < 2. The systematic uncertainties associated with the combined results are fully correlated bin- 
to-bin in each distribution. They are approximately 3% at midrapidity, and rise to about 10% at forward 
and backward rapidity. 


4 Results 

4.1 Z -* tt cross section 

From the combined Z —> ee and Z —> pp data a total cross section of 139.8 ± 4.8 (statistical) ± 6.2 
(systematic) ± 3.8 (luminosity) nb is obtained in the \y* z \ < 3.5 acceptance. Based on the MC simulation 
(and the models discussed below) this acceptance covers approximately 99.5% of the total Z —> €£ cross 
section. Restricting the results to the smaller rapidity interval of -3 < y* z < 2, the cross section is 119.3 ± 

2.2 (stat.) ± 3.4 (syst.) ± 3.2 (lumi.) nb. Table 3 lists the integrated cross section in the larger and smaller 
rapidity ranges as measured for each channel and their - combination. 


* 

y z 

[-2,0] 

[0,2] 

[-3,2] 

[-3.5,3.5] 

Z-^pp 

54.2 ± 1.6 ± 1.3 

45.3 ±2.1 ±0.9 

118.2 ±3.3 ±2.6 

N/A 

Z —> ee 

55.1 ± 1.8 ±5.9 

46.5 ±2.2 ±5.0 

121 ±3 ± 13 

143 ±5 ± 17 

Z->M 

54.4 ± 1.3 ± 1.4 

45.9 ± 1.4 ± 1.4 

119.3 ±2.2 ±3.4 

139.8 ±4.8 ±6.2 

CT10 (NLO) 

47.4 ± 0.9 

46.8 ± 0.9 

110.8 ±2.9 

132.2 ±3.3 

CT10+EPS09 (NLO) 

48.7 ± 1.0 

43.5 ± 1.1 

108.6 ±3.1 

127.4 ± 3.6 

MSTW2008 (NNLO) 

48 3 +1-2 

-0.9 

47 9 +1-2 

+ '-0.9 

113.5+11 

135.2^2 7 


Table 3: The measured integrated cross section, in nb, for several rapidity ranges, for Z —* pp, Z —> ee, and the 
combined Z —» CC. The first uncertainty listed is statistical and the second systematic. There is an additional 2.7% 
luminosity uncertainty on each cross section. The cross sections predicted by the models (see text) are also shown. 
The uncertainties listed with the model calculations are the PDF and scale uncertainties added in quadrature. 


The measured cross section may be compared to a p+ Pb model prediction composed of a linear sum of the 
nucleon-nucleon cross sections: 82 cr{pp —> Z+X ) + 126 cr(pn —> Z+X), corresponding to the numbers of 
protons and neutrons in the Pb ion. The value of cr(pn —> Z + X) is 2% higher than that of cr(pp —> Z + X) 
in all models discussed below. Calculating the baseline nucleon-nucleon cross sections using the CT10 
PDF at next-to-leading order (NLO), as in the corresponding MC simulation, the model yields values of 

132.2 ± 3.3 nb in the range \y z \ < 3.5, and 110.8 ± 2.9 nb for -3 < y* z < 2 where the uncertainties 
are the sums in quadrature of PDF and scale (renormalization and factorization) uncertainties. Using the 
MSTW2008 PDF, calculated with Fewz [37] at next-to-next-to-leading order (NNLO), cross sections of 

135.2 nb are obtained for \y* z \ < 3.5 and 113.5 ^'2 f° r - 3 < y z < 2. At NLO the results from 
MSTW2008 are very close to the CT10 results. In addition to the simple model of the /;+Pb Z boson 
cross section as a linear sum of nucleon-nucleon cross sections, calculations are performed incorporating 
nuclear corrections of the PDF. Including the EPS09 modifications [38] to the CT10 PDF results in cross 
sections of 127.4 ± 3.6 nb and 108.6 ±3.1 nb, respectively. 

For a more detailed understanding of Z boson production, the measured cross section as a function of the 
Z boson rapidity is presented in Figure 4 and compared to model calculations. The data are seen to be 
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strongly asymmetric about y* z =0. The CT10+EPS09 calculations come closest to reproducing the shape 
of the measured y* z differential cross section. A x 2 test of compatibility between the data and the model 
shapes (irrespective of normalization) finds that the CT10+EPS09 shape of the y* z distribution gives a 
p-value of 0.79. The unmodified CT10 calculation and MSTW2008 calculations have /^-values of 0.07 
and 0.01, respectively. A Kolmogorov-Smirnov test was also performed and resulted in probabilities of 
0.96, 0.09, and 0.07 for CT10+EPS09, CT10, and MSTW2008 model calculations. This is consistent 
with the preference for the observation of nuclear correction effects as in the x 2 test. 



_J_._i_._i_._._._i_ l~ 

-2 0 2 y* z 


Figure 4: (Color online) (a) The dcrldy* z distribution from Z — * ft, shown along with several model calculations in 
the upper panel. The bars indicate statistical uncertainty and the shaded boxes systematic uncertainty on the data; 
the uncertainties on the model calculations are not shown, (b-d) Ratios of the data to the models. The uncertainties 
of the model calculations (scale and PDF uncertainties added in quadrature) are shown as bands around unity in 
each panel. An additional 2.7% luminosity uncertainty on the cross section is not shown. 


Nuclear modification of PDFs is fundamentally related to the Bjorken x of the relevant parton. At leading 
order, x p in the proton and xpb in the lead nucleus are related to the reconstructed Z boson kinematics by: 

m a e y z m u e~ lJ z 

x p = - , x Pb = —-= . (1) 

V^NN V^NN 

The resulting xpb distribution is shown in Figure 5 and compared to model calculations. 

Figure 6 shows the distributions for -3 < y* z < 2, and separately for -2 < y* ? < 0 and 0 < y* 7 < 2. 
These are compared to the baseline CT10 model. The dependence is less sensitive to nuclear effects 
and good agreement between the experimental measurement and the MC simulation shape is observed. 
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Figure 5: (Color online) (a) The differential cross section of Z boson production multiplied by the Bjorken x of the 
parton in the lead nucleus, xpbdcr/dxpb, as a function of xp h using Z —> (( events shown along with several model 
calculations in the upper panel. The bars indicate statistical uncertainty and the shaded boxes systematic uncertainty 
on the data; the uncertainties on the model calculations are not shown, (b-d) Ratios of the data to the models. The 
uncertainties of the model calculations are shown as bands around unity in each panel. There is an additional 2.7% 
luminosity uncertainty on the cross section. 


4.2 Centrality-dependent yield 

Results are presented for the centrality-dependent Z boson yield. If the rate of Z boson production were 
consistent with geometric expectations, then the Z boson yield divided by (A co ii> should be independent 
of centrality. To investigate this, the yield of Z bosons per event scaled by {N co \\), within -3 < y* z < 2, 
is displayed as a function of (/V part ) in Figure 7. The yield is independent of centrality defined using 
the standard Glauber model. Using the GGCF centrality models increases (N co n) in central events and 
reduces it in peripheral events, consequently the yield divided by (iVcoii) is reduced in central events and 
increased in peripheral events. Figure 7 also shows the yield without the application of the centrality bias 
corrections discussed in Sec. 3.3. 

The ATLAS Collaboration has previously measured the inclusive charged-hadron multiplicity in p +Pb 
collisions as a function of centrality [22], and the centrality dependence of that quantity is similar to that 
observed in the present measurement. In order to quantify the similarity, the ratio (dNz / dy* z ) /(<7(V c h/ drj) is 
plotted vs (Apart) in Figure 8. The charged-particle yield is expected to scale with (/V parL ) and the Z boson 
yield with (A co ii) = (/V par[ ) - 1, and so the ratio is fit to a function with the form a ■ ((A part ) - l)/(A part ). 
This function describes the data well for the GGCF cases, and less so for the standard Glauber model. 

To further investigate the behavior observed in the rapidity differential cross section, the y* z dependence 
of the Z boson yield in different centrality bins is also measured, as shown in Figure 9. The differences 
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1 10 10 2 pz [GeV] 1 10 10 2 [GeV] 

Figure 6: (Color online) (a, c) The distributions of the differential cross section of Z boson production as a function 
of the transverse momentum of the Z, shown along with the CT10 model calculation, (b, d) Ratio of the data to 
the model. (a,b) for -3 < y* z < 2; (c,d) for -2 < y* z < 0 and 0 < y* z < 2. The bars indicate statistical uncertainty 
and the shaded boxes systematic uncertainty. The leftmost bin represents the range 0-3 GeV. An additional 2.7% 
luminosity uncertainty of the cross section is not shown. 


between data and the model are larger in central collisions. The (_/V co u)-sealed ratio of central to peripheral 
data, Rqp, defined as 

n , *, <iVcoii) peripheral dN z ential /dy* z 

Rcp(y z ) = ——:—— x 


<AUi>' 


central 


dm 


peripheral 


( 2 ) 


•Z /dy* z 

is used to observe changes in the rapiditiy distribution for different centrality bins in a model-independent 
way and is shown in Figure 9. Events with 40-90% centrality define the peripheral event selection, and 
two central selections, 0-10% and 10-40%, are compared with it. A linear fit of the Rqp (y* z ) for 0-10% 
centrality, results in a slope of -0.11 ± 0.04, which suggests that the y* z distribution may be different in the 
most central events compared to peripheral events. For 10-40% centrality, the slope is -0.05 ± 0.03. 
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Figure 7: (Color online) The yield of Z bosons per event scaled by the mean number of nucleon-nucleon collisions 
(iVcoii) as a function of the mean number of participating nucleons (N part ). Each panel uses a different Glauber 
model configuration in calculating ( N palt ) (and (/V co n)). (a) Standard Glauber model with no Glauber-Gribov color 
fluctuations; (b) GGCF with = 0.11; (c) GGCF with o.) rr = 0.2. The data are compared to the CT10 model 
prediction plotted at ( N palt ) = 2. The bars indicate statistical uncertainty and the shaded boxes systematic uncer¬ 
tainty. The systematic uncertainties are correlated bin-by-bin. The (/V co n) uncertainty plotted does not include the 
bin-by-bin fully correlated uncertainty stemming from the uncertainty on cr(N + N —» X), which is instead included 
in the CT10 prediction uncertainty. As a reference, the data are plotted as they would be with no centrality bias 
correction in the open points. 
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Figure 8: (Color online) The ratio of the Z boson multiplicity to the inclusive charged particle multiplicity, 
(dNz / dy* z ) / (dN c h/chj), as a function of (fV part ). A function of the form a ■ ((fV part ) - l)/(A^ palt ) is also shown. The 
normalization a is set based on the GGCF with Wo- = 0.11 points. Statistical uncertainties are plotted as bars and 
systematic uncertainties as shaded boxes. 



Figure 9: (Color online) (a) The rapidity differential Z boson yields, scaled by (/V co ii), for three centrality ranges 
compared with the CT10 model calculation. The bars indicate statistical uncertainty and the shaded boxes system¬ 
atic uncertainty. The (/V co n) is defined with the standard Glauber model (ov=0). The scale uncertainty stemming 
from the centrality calculation for each bin is included in the systematic uncertainty. The uncertainty associated 
with the model is not plotted, (b) The ratios of the data to the model. The uncertainty of the model added in quadra¬ 
ture with the scale uncertainty due to uncertainty in the inclusive NN cross section is shown as a band around unity, 
(c) Rep (see text for details). The 0-10% and 40-90% centrality points are offset for visual clarity. The arrows in 
the lower panels indicate values outside the plotted axes. 
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5 Summary 


The Z boson production cross section has been measured in p +Pb collisions at yJs NN = 5.02 TeV with the 
ATLAS detector at the LHC, using Z —» ee (Z —» pp) decays in a 29.4 ntr 1 (28.1 nb 1 ) data sample. It 
is found to be slightly higher than predictions based on perturbative QCD calculations. Disregarding the 
difference in overall normalization, the shapes of the y* z and .rpt,- dependent cross sections are somewhat 
better described by models that include nuclear modification of the lead nucleus PDF compared to those 
that do not, although models without nuclear modification are not excluded. Following the application of 
a centrality bias correction, the centrality-dependent yield is found to scale with (N c0 \\). In addition, the 
centrality dependence of the y* z distribution was studied, and the asymmetry in y* 7 was found to be slightly 
larger in more central events. Integrated over y* z , the centrality dependence appears to be consistent with 
binary scaling and is similar to the production of inclusive charged particles. 
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